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Abstract. We investigate the spectral distribution of the damped wave equation 
on a compact Riemannian manifold, especially in the case of a metric of negative 
curvature, for which the geodesic flow is Anosov. The main application is to 
obtain conditions (in terms of the geodesic flow on X and the damping function) 
for which the energy of the waves decays exponentially fast, at least for smooth 
enough initial data. We review various estimates for the high frequency spectrum 
, in terms of dynamically defined quantities, like the value distribution of the time- 

(N ■ averaged damping. We also present a new condition for a spectral gap, depending 

^ |. on the set of minimally damped trajectories. 

D 

1. Introduction 

(-1 

p_il 1.1. Spectrum of the damped wave equation. Given a Riemannian manifold 

{X, g) and a damping function a G C°°{X, M+), we are interested in the solutions of 
the damped wave equation (DWE) 



1.1) {Of - A + 2a{x)dt)v{x, t) = 0, v{x, 0) = vq, dtv{x, 0) = vi . 



The natural setting is to take intial data (fo,fi) in the space Ti *== H^{X) x L'^{X) 
^ ■ This equation is equivalent with the system 

(1-2) {^^, + A)v{t) = 0, ^=(_°^ _2^^), viO)^^ 

o , 

! with the correspondence v{t) = {v{t),idv{t)). A generates a strongly continuous 

O I semigroup on "H, so the solution to fll. 1111.21) reads 

^ (1.3) v{t) = e-''^v{0). 

We will always assume that the damping is nontrivial, a ^ 0. Apart from the 
^ ■ constant solution v{t) = (1, 0), all solutions then decay. Physically, this is expressed 

by the fact that the energy of the waves, 

(1-4) E{v{t)) = l{\\Vv{t)f + \\dMt)f), 

will decay to zero when t — )■ oo for any v{0) G Ti. In some sense, the waves are 
stabilized by the damping. 

To analyze this decay, it it natural to try to expand the solution in terms of the 
spectrum of A. This spectrum is discrete, consisting of countably many complex 
eigenvalues {r„} with Rer„ — )■ ±oo. It can be obtained by solving the generalized 
eigenvalue equation 

(1.5) (-A - - 2zar) u = 0. 
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Figure 1.1. A damped geodesic. 



The following properties are easily shown {a^in = miux^zx a(x), similarly for amax)- 

Proposition 1. jLeb93j 74// eigenvalues except tq = satisfy Imr„ < 0. 
// Re r„ 7^ then - Im r„ G [a^in, Omax] ■ 
The spectrum is symmetric w.r.to the imaginary axis. 

To each eigenvalue r„ corresponds a quasi-stationary mode Un{x), an eigenstate 
Un = {un, TnUn) of A rcads, and a solution Vnit, x) = e~**'^"M„(x) of the DWE. Hence 
Imr„ represents the quantum decay rate of the mode m„. 

If flmin > 0, then the energy decays exponentially, uniformly for initial data v{0) G 
"H. Precisely, there exists 7 > 0, C > such that 

(1.6) E{v{t)) <Ce-^^' E{v{0)), Vv(0) G Vt > 0. 

1.2. The Geometric Control Condition. The condition amin > is not neces- 
sary to ensure such a uniform exponential decay. Since all eigenvalues except tq = 
satisfy Imr„ < 0, for any C > the subspace Tic C "H spanned by the eigenstates 
{un, |Rer„| < C} is finite dimensional, and for any initial data v{0) C He the en- 
ergy will decay exponentially. Hence, the failure of exponential decay can only come 
from the behaviour of waves at high frequency. In this high frequency limit, a natu- 
ral connection can be made with the classical ray dynamics on X, equivalently the 
geodesic flow $* on the unit cotangent bundle S*X (see §1.4p . Using this connection. 
Ranch and Taylor |RauTay75] showed that the uniform exponential decay fll.6p is 
equivalent with the Geometric Control Condition (GCC), which states that every 
geodesic meets the damping region {x G X, a{x) > 0} (due to the compactness of 
S*X each geodesic does it within some time Tq > 0). 

This condition can be expressed in terms of the time averages of the damping, 
namely the functions 

1 /■* 

(1.7) {a)t{p) = - ao$^(p)ds, pGT*X, t>0. 

Jo 

GCC is equivalent to the fact that, for t > large enough (say, t > 2Tq), the function 
{a)t is strictly positive on S*X. Lebeau |Leb93j generalized this result to the case of 
manifolds with boundaries. He also showed that the optimal decay rate 7 is given 
by min (G, a_), where G = inf {— Imr„, r„ 7^ 0} is the spectral gap, while 

(1.8) a_ == lim min{a)t 

t^oo S*X 

is the minimal asymptotic damping. We will explain the relevance of a_ in §2.2[ 

Koch and Tataru [KoTa94j studied the same question in a more general context 
(case of manifolds with boundaries, and of a damping taking place both in the 
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"bulk" and on the boundary). They showed that the averages {a)t govern the decay 
of the semigroup, up to a compact subspace. 

Theorem 2. |KoTa94| Thm 2] For each e > and each t > there exists a subspace 
'He,t C l-L of finite codimension such that 

< exp |-tmin(a)i} + e. 

For any subspace "Hi C "H of finite codimension, 

exp |— t min(a)t| < 



-itA\ 



A consequence of this result is the characterization of the Fredholm spectrum of 
the semigroup. In the present situation of damped waves on a compact X withough 
boundary, their result state£| that this Fredholm spectrum is given by the annulus 
{z e C, e-*"+ < 1^1 < e"*"-}. 

Notice that it is possible to have a positive gap G > and at the same time a_ = 
(failure of the GCC), see e.g. an example in |Ren94j . This reflects the fact that the 
spectrum of the semigroup e~**'^ is not controlled by the spectrum of generator A, 
a frequent problem for nonnormal generators like A. 

1.3. Beyond Geometric control: a few cases of uniform energy decay for 
regular data. If GCC fails, one can construct initial data v{0) G H such that 
E{v{t)) decays arbitrarily slowly |Leb93] . Yet, it is possible to show that more 
regular data, say in some Sobolev space T-L'^ = H^~^^ x if"*, s > 0, decay in some 
uniform way. In the most general situation, Lebeau showed a logarithmic decay of 
the energy |Leb93t Thm 1], 

E(vm" < g. '°^'f + '°!'V'" Mn. . V. = .(0) e n; w > o, 

log(3 + 1)* 

The proof uses Carleman estimates, which imply some control of the resolvent 
{t — A)^'^ in a strip {Rer > C, |Imr| < /(Rer)} for some function /(r) which 
decays exponentially as r — )■ oo. By interpolation, proving such an estimate for 
some So > implies a similar estimate for all s > 0. 

In specific dynamical situations, one can show resolvent estimates in a larger strip, 
e.g. for some function /(r) decaying algebraically: one then gets an algebraic decay 
0{t~'^'' \\v\\^s) for the energy of regular data. Burq and Hitrik showed that this is 
the case for damped waves in the (chaotic) stadium billiard, if the damping function 
vanishes only in some part of the central rectangle, so that the set of undamped 
trajectories, 

(1.9) K = {p G S*X, Vt G M, a($*(p)) = O} , 

consists in a collection of (neutrally stable) "bouncing ball" trajectories jBuHi07j . 

Christianson studied the case where the set K consists in a single hyperbolic closed 
geodesic: the energy decay is bounded by a stretched exponential 0{e^^'''^^^'^ \\v\\y^s) 
|Chris07[ Thm 5'] [ChrisTT] . The hyperbolicity of the geodesic induces a strong 
dispersion of the waves, responsible for this fast decay. 



under the condition that, for any T > 0, the set of T-periodic geodesies has measure zero. 
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In §3.31 and 13.41 we will present situations, with {X, g) a manifold of negative 
curvature, for which the energy of regular data decays exponentially: 

(1.10) E{v{t)y/'^ < Cs e-^^' \\v\\^, , G Vt > 0. 

The proof proceeds by controlling the resolvent (r — A)~^ in a strip {Re t > C, | Im r | < c}, 
and then uses standard arguments |Chris09j . Here as well, the resolvent bounds are 
based on certain hyperbolic dispersion estimates. See the Corollaries [16] and [19] for 
the precise dynamical conditions, which involve the interplay between the flow and 
the damping function. 

1.4. High frequency limit — semiclassical formulation and generalization. 

As explained above, the decay of the energy is mainly governed by the high fre- 
quencies Rer ^ 1. A semiclassical formulation of the problem was used in |SjoOO| , 
allowing to use /i-pseudodifferential techniques. Take an effective Planck's "con- 
stant" < ^ <^ 1. In order to study the frequencies Rer = + 0{1), we perform 
the rescaling 

(1.11) r = ^ with z in the disk D{l/2, CK). 
The equation (11.51) becomes 



[1.12) {P{h,z)~z)u = 0, P{h,z) = -'-^-ihV2za = -'^ iha + 0{h^). 

More generally, we consider operators of the type 

(1.13) p(n,z) = -^ + thOprU)- 

where the symbofl G S°{T*X) is now a function on the phase space T*X, and 
depends holomorphically on z ^ D{l/2, Ch), with q^ = qi/2 + 0{h) |SjoOO| . With a 
slight abuse, will call q = 51/2(2;, "the damping function" (a damping effectively 
occurs only at phase space points where q{x,C,) < 0, while positive values of q rather 
correspond to a "gain" or "pumping"). 

In this framework, we are interested in the semiclassical distribution of the gen- 
eralized eigenvalues Zn{h) of the operator P{h,z) in the disk D{l/2,Ch). Because 
we will specify this spectrum at the scale we may for convenience drop the term 
0{h^) in (11.131) . and rather study the spectrum of the z-independent operator 



:i.l4) P{h)=Po{h)+ihOvn{q). Po{h) = 



The principal symbol of this operator is Po{x,C, ) = it generates the geodesic flow 
$* = exp(tHpg) on T*X, with unit speed on the energy shell Pq^{1/2) = S*X. 



^For any A: G K, the symbol class S^{T*X) denotes the functions a(a;, ^; h) satisfying the follow- 
ing bounds; for any multi-indices a,f3 £ N'', 

d^d^aix, h) < C^^^ (1 + |C|)'=-I''' uniformly for h G (0, 1), 

By Op;j((7) we refer to a semiclassical quantization on X, with the property that Op;j((j') is sym- 
metric if q is real (see e.g [EvZwi App. E]). The operators 0\>ySsI) form the class "^^(X) of 
pseudodifFerential operators (PDO). 
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2. Semiclassical spectral distribution of P{h). General results 



The results we present in this section hold for arbitrary Riemannian manifolds 
{X,g). We skip the proofs. 

2.1. Weyl law for the real parts. In the semiclassical limit, the number of 
eigenvalues of the selfadjoint operator Po{h) = — in the interval [| ± e] == 
[| — e, I + e] is approximately given by Weyl's law. When applying the perturba- 
tion ihOpfi{q) = 0{h), it seems natural to expect that most eigenvalues of Po{h) 
in this interval will remain inside a (complex) neighbourhood [| ± e] + 0{h). One 
can indeed show that the "macroscopic" spectral distribution is given by the same 
Weyl's law as in the undamped case. 

Theorem 3. |MarMat84l |SjoOO] For any e > small enough (possibly depening on 
h), the number of eigenvalues of P{h) in the strip = {| — e < Re^; < | + e} is 
semiclassically given by 



# {SpecP(n) n S,} = i2nh)-'^Vol {p E T*X, \po{p) - 1/2| < e} + O (r'^+i) . 

Here Vol corresponds to the symplectic volume in T*X . This estimate implies that, 
for C > large enough, 

(2.1) # {Spec P(n) nD(l/2, C;i)} X h-^+\ 

This result is actually not so easy to prove. It uses some subtle analytic Fredholm 
theory in order to relate the spectrum of P{h) with that of a simpler "comparison 
operator" . 

2.2. Restrictions for the imaginary parts (quantum decay rates). Even- 
though the operator P{h) was inspired by the wave equation, it is convenient for us 
to consider it instead as the generator of a Schrodinger equation 

(2.2) ihdtv - P{h)v = 0. 

Any eigenvalue/ vector (z„,m„) of P{h) is thus associated with a solution Vn{x,t) = 
e~*^"*/^ti„(x) of this equation, and we will also call Im Zn/h the quantum decay rat^ 
associated with m„. 

First we remark that an eigenmode m„ associated to an eigenvalue Zn € D{l/2,Ch) 
is semiclassically microlocalized on ^(1/2) = S*X. To state this property we need 
to introduce energy cutoffs. 

Definition 4. An energy cutoff will be function x ^ ([1/2 ± 5] , [0, 1]) for some 
(5 > 0, and satisfies x{^) = 1 iisar s = 1/2. The cutoff Xo is said to be embedded in 
Xi, and we note xi ^ Xo, iff Xi = 1 ii^ar suppxo- From the cutoff x we construct 
the (pseudodifferential) cutoff operator x(Po(^)) G \E'~°°(X). 

Proposition 5. Take an energy cutoff x- Let u = u{h) be an eigenstate of P{h) 
with eigenvalue z{K) G D{l/2,Ch). Then, 

\\ii~xiPo))u\\ = Oih°^)\\u\\. 

This localization property explains why only the neighbourhood of S*X will be 
important for us. 

"'From (|l.lip we have Im z„/fi, = ImT„+C'(?i), so this denomination is (approximately) consistent 
with the one of m.l[ 

5 



2.2.1. A factorization of the Schrddinger propagator. An easy way to realize how 
the skew-adjoint subprincipal term ihOpfj{q) imphes a "damping" is to analyze the 
propagator for the Schrodinger equation (12. 2p . The following Proposition adapts a 
more general result due to Ranch and Taylor |RauTay75] . 

Proposition 6. Assume q G C'^{T*X), and take P{K) as in 1^1. 1^^ . For any fixed 
t G M, decompose the Schrddinger propagator = e"**-^/^^ into 

(2.3) = f/* B{t) , where f/* = e"'*^"/'* is the undamped propagator. 

The operator B[t) is a PDO in \&°(X) of principal symbol 

hit, p) = exp ( / go $^(p) ds\ = exp {t {q)t{p)} , 



where the time averaged damping {q)t generalizes ( [i. 7| j. 

The factor b{t,p) is the accumulated damping along the orbit {$'^(p), < s <t}. 
If one starts from a Gaussian wavepacket uq microlocalized at a point p, then the 
state V^uq is a wavepacket microlocalized at $*(p), and the above statement shows 
that the norms are related through this factor: 

\\ut\\ = \\uo\\ {b{t, p) + 0{h)) . 

Applying an energy cutoff allows to extend the factorization (12.31) damping functions 
q G S^(T*X) of noncompact support. 

Lemma 7. Consider P{h) with a damping function q G S^{T*X). Take two embed- 
ded cutoffs Xi Xo supported in [1/2 ± 6], and the truncated damping q =^ Xi{Po) <?■ 
Then, for any t > fixed, one has 

(2.4) V'xoiPo) = V'xoiPo) + Ol^^l^K^) , 
where is the propagator corresponding to P = Pq + ihOp^{q). 

This identity uses the fact that the propagation does not modify the energy 
localization properties. 

2.2.2. From propagator factorization to a resolvent bound. We can now easily obtain 
a first constraint on the quantum decay rates. Applying Proposition[6]to the function 
q = Xi{Po)(l of Lemma [71 we get for any fixed t > 0: 

\\V'xoiPo)\\^.^L.= Bit)xoiPo) ^^ + Oih^) 

(2.5) = maxbit)xoiPo) + Otih) 



< exp < t max {q)t > + Ot{h). 

{ suppxoopo J 

On the second line we used the sharp Garding (in)equality for the operator B{t) G 
\E'°(X). The maximum maXp-i^j^^2±5]) (Q')t decreases when t — ?■ oo or when 5 — )■ 0, 
and converges to the asymptotic maximum on S*X, 

(2.6) q+=\immax{q)t{p). 

t^oo peS*X 

One similarly defines an asymptotic minimum q_. We then get the following norm 
bounds for the propagator. 
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Proposition 8. Fix e > 0. If the energy cutoff xo has a small enough support, and 
> is large enough, then for any (fixed) t>T^ and any small enough h> 0, the 
following bounds holds: 

(2.7) \\V'xo{Po)\\l.^l. < e^'^+'^\ \\V-'xo{Po)\\l.^l. < e^-^-^^^ 
From there one easily obtains the following resolvent and spectral bounds. 

Theorem 9. |Leb93t SjoOO| Take any e > 0. Then, there exists h^, > such that 
for h < h^, the following resolvent estimate holds: 

(2.8) \/z E D{l/2,Ch)\{z, Imz/he [Q- - e, q+ + e]} , \\iPih) - zy^\\ < ^. 

As a consequence, for h < all eigenvalues z„ G D{l/2,Ch) of P{h) satisfy 

Im z„ 



- G [g_ - e, 



h 

In the case of the damped wave equation {q{x,S,) = —a{x)) and in a situation of 
geometric control (a_ = — g+ > 0), the resolvent estimate (12. 8p can be used to show 
the uniform exponential energy decay (II. 6p |Hit03j . 

2.3. Questions on the spectral distribution. So far we showed that the quan- 
tum decay rates are bounded by the asymptotic extrema of the time-averaged damp- 
ing. The following questions were raised in |SjoOO lAschLebj lAnanlOj concerning 
their semiclassical distribution. 

(1) What are the possible accumulation points of the quantum decay rates when 
h ^ 07 In particular, are there sequences of decay rates {lmz{h)/h)fi^Q 
converging to the extremal values q±7 

(2) Do the quantum decay rates admit an asymptotic distribution when 0? 
Namely, for a given interval / C [q-, q+] and 1 ^ e{h) ^ h, does the ratio 



G SpecP(;i), 


Rez- 1/2| < e 


Im z ^ t\ 


#{2 G SpecP(;i), 


Rez - 1/2 





have a limit when h ^ 07 Is this limit distribution related with the value 
distributions of the averages 



3. Spectral estimates on Anosov manifolds 

The above questions are open in general. In order to get more precise informations 
on the spectrum of P{h), one needs to make specific assumptions on the geodesic 
flow on X. For instance, the case of a completely integrable dynamics has been con- 
sidered by Hitrik-Sjostrand in a sequence of papers (see e.g. |HitSjo08] and reference 
therein). The case of nearly-integrable dynamics including KAM invariant tori has 
been studied by Hitrik-Sjostrand-Vii Ngoc jHSVNOT] . In these cases, one can trans- 
form the Hamiltonian flow into a normal form near each invariant torus, which leads 
to a precise description of the spectrum "generated" by this torus. A Weyl law for 
the quantum decay rates was recently obtained in [Hit Sj oil] (for skew-adjoint per- 
turbations i9{h) Op;j(g), with 9{h) ^ h). On the other hand, Asch-Lebeau [AschLeb] 
addressed Question 1 for the case of the 2-dimensional standard sphere. They show 
that, if the damping function q has real analytic real and imaginary parts, there is 
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Figure 3.1. Structure of the Anosov flow near an orbit $*(p) 



generically a spectral gap 7 > 0: for h small enough, all eigenvalues Zn G D{l/2,Ch) 
have quantum decay rates 

Im Zn/ he[q-+ 7, g+ - 7] • 

The proof proceeds by applying a complex canonical transformation, such that the 
real part of the pulled-back damping function takes values in the above interval. In 
this case, the range of the quantum decay rates is strictly smaller than the range of 
classical decay rates. 

We will see below that such a spectral gap may also occur in the case of Anosov 
geodesic flows, thanks to a different mechanism, namely a hyperbolic dispersion 
property due to the instability of the classical flow (see Thms fTSl and fTSl below) . 

In the next section we recall the definition and properties of Anosov manifolds. 

3.1. A short reminder on Anosov manifolds |KatHas95] . At the "antipode" 
of the completely integrable case, one finds the "strongly chaotic" flows, namely the 
Anosov (or uniformly hyperbolic) flows. Uniform hyperbolicity means that at each 
point p G S*X there exists a splitting of the tangent space, 

TpS*X = RHp{p) © E+{p) © E-{p), 

where Hp{p) is the Hamiltonian vector field, E'^{p) are the unstable and stable 
subspaces at the point p; they both have dimension (i — 1, and are uniformly trans- 
verse to e.o.. The families {E^{p), p G S*X} form the unstable/stable distributions, 
they are invariant w.r.to the flow. Holder continuous, and are characterized by the 
following property: there exists C, A > such that 

(3.1) Vp e S*X, Vt; G E^(p), Vt > 0, \\d^fv\\ < C e~^' \\v\\ . 

As was shown by Hadamard, such a geodesic flow is obtained if the manifold (X, g) 
has a negative sectional curvature. The long time properties of such a strongly 
chaotic flow are well understood. In particular, Hopf proved that this geodesic flow 
is ergodic w.r.to the Liouville measure on S*X. 

Definition 10. Call pi the normalized Liouville measure on S*X. Then, the ge- 
odesic flow on S*X is ergodic (w.r.to pi) iff, for any continuous function q on 

S*X, 

for Pl — almost every p G S*X, lim {q)t{p) = q'= Piiq)- 
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The following quantities provide quantitative measures of the hyperbolicity: 



(3.2) the maximal expansion rate Amax *== lim sup - log \ \d^ 



the unstable Jacobian J^(p,t) = det \e+{p)) , t > 0, 

(3.3) and its infinitesimal version Lp~^{p) —log J^{p,t) \t=o, 

(3.4) the minimal expansion rate i^min'= lini min - log J^(p, t) . 

a — 1 t^oo pes*x t 

The Jacobian depends on the choice of norms on the spaces TpS*X. However, its 
long time asymptotics is independent of this choice. In variable curvature these 
rates are related by Amax > t'min > A > 0. In particular, if the positive Lyapunov 
exponents are not all equal, one has v-rain < Amax- 

A particular class of Anosov manifolds consists in quotients of the dimensional 
hyperbolic space by co-compact subgroups of its isometry group. These manifolds 
have a constant curvatur^ — A^. Using the natural norms on TpS*X, one then has 

J+(p,t) = e*('^-l)'^, y.+ (p) = A(rf- 1), Amax = J^min = A. 

In this case the contraction/expansion are both homogeneous (independent of the 
point p) and isotropic (independent in the direction). 

Let us finally notice that the study of the operator P{h) on an Anosov manifold 
belongs to the field of "quantum chaos". The methods we will use below occur in 
various problems of this field, e.g. the study of eigenstates of the Laplacian Po{h) 
on such manifolds |Zel09] . 

3.2. Fractal Weyl upper bounds for the quantum decay rates. In this section 
we address Question 2, that is the asymptotic distribution of the quantum decay 
rates, for the case of Anosov manifolds. 

3.2.1. Typical quantum decay rates for ergodic flows. A basic property of the geo- 
desic flow on any Anosov manifold is ergodicity. Sjostrand has studied the asymp- 
totic disribution of the quantum decay rates for any ergodic geodesic flow |SjoOO| , 
and obtained the following result. 

Theorem 11. |SjoOO| Assume the geodesic flow on S*X is ergodic w.r.to the Liou- 
ville measure. Then, For any C > 0, e > 0, one has as h ^ 0, 

# |z G SpecP(n) n 1^(1/2, Ch), ^ [g - e, g + e]| = o{h-''+'). 

Comparing this bound with the Weyl law (12. ip shows that, in the semi classical 
limit, almost all the quantum decay rates are close to the phase space average g. 
The latter could be called the typical value for quantum decay rates, while values at 
finite distance from g are atypical. The asymptotic distribution is simply the delta 
measure 6q. 

Remark 12. As was noticed in [AnanlOj . this concentration result could be seen 
as a "nonperturbative version" of quantum ergodicity. The latter property states 



^Usually one normalizes the metrics on H'^ so that the curvature is —1. We prefer to keep track 
of the curvature in our notations. 
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that almost all the eigenstates m° of Po(^) with eigenvalues z° G D{l/2,Ch) sat- 
isfy {Un,Op^{q)u'^) = q + on^oi^l). A naive perturbation theory argument would 
predict that, when switching on i^Op^j(g), the eigenvalues z° move to Zn = + 
ih{u^, Op;j(g)M°)+C(/i^), so that the quantum decay rates Ivuzn/h = Op^{q)u'^) + 
0{h), the RHS being equal to g + o(l) for almost all n. Of course, this argument 
does not apply because the perturbation z/iOp^j(g) is much stronger than the mean 
spacing between successive eigenvalues. 

To prove the above theorem one relates the counting of quantum decay rates with 
the value distribution of the quantum averages {q)t on S*X, that is the volumes 

(3.5) HL{pe S*X, {q)tip)>a}, a G M, t > 0. 

Indeed, the main intermediate result in the proof is the bound 
(3.6) 

{I in z 1 
zeSpecP{h)nD{l/2,Ch), > a \ < Ct,5h-''^' Pl{p e S*X, {q)t{p)>a- 

which holds for any fixed t > and 6 > (this bound holds independently of the 
ergodicity assumption). Ergodicity then implies that the value distribution of {q)t 
converges to 6g when t — >■ oo; in particular, for any a > q the volume (13. 5p decays 
to zero when t ^ oo. We have then obtained the bound o{h~'^~^^) for the quantum 
decay rates > a > q. The case of values < a < g is treated analogously. □ 

3.2.2. Large deviation estimates for Anosov flows. For an Anosov flow, one has more 
precise estimates for the volumes (13. 5p . which can then induce sharper bounds on 
the number of the atypical quantum decay rates. These volume estimate take the 
form of large deviation estimates. Let us introduce some notations. We call Ai 
the set of $*-invariant probability measures on S*X. For each such measure p, we 
denote by h^sip) its Kolmogorov-Sinai entropy: this is a nonnegative number, which 
measures the complexity of a /i- typical trajectory |KatHas95] . Then, we define the 
rate function : M — )■ M as follows: 

(3.7) VsGM, H{s) = ^nv{hKs{p)-p{^'^),peM,p{q) = s], 

where we recall that is the infinitesimal unstable Jacobian (13. 3p . We are now 
ready to state our large deviation result. 

Theorem 13. IKifQOj ylggnme the geodesic flow on S*X is Anosov. Then, for any 
closed interval / C M and for any q G C°°{S*X), the time averages {q)t satisfy 

(3.8) limsup^log/iL{pG {q)t{p) e 1} < sup H{s). 

The rate function H is continuous on [g_,g+], smooth and strictly concave on 
(g_,g+), negative except at the point H{q) = 0, satisfies H{s) > — sup ^^j^ p{ip~^) 
for s G [g_,g+], and is equal to — oo outside [g_,g+]. As a consequence of (13. Sp . for 
any a > q and e > arbitrary small, there exists Ta^e > such that 

(3.9) Vt > T,,e, PL {p e S*X, {q)t{p) >a}< e*(^(°)+^). 

Due to the negativity of H{a) for a q, we see that the probability of {q)t taking 
atypical values decays exponentially when t — )■ oo. 
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3.2.3. Fractal Weyl upper hounds on Anosov manifolds. Using these large deviation 
estimates, Anantharaman [AnanlOj improved Thm[TT]by letting the averaging time 
grow with h'm & controlled way. The optimal time is the Ehrenfest time 

(3.10) T = T^,. = (l-2e)^°^^/^ 



A 



max 



where A^ax is the largest expansion rate f l3.2p and e > arbitrary small. What is the 
signification of this time? For any / G S~°°{T*X) supported in an e-neighbourhood 
of S*X, the classically evolved observable /o$* remains in the "good" symbol clas^ 
Sy^_^{T*X) uniformly for times |t| < TEhrf^- In turn, the symmetric averages 

(3.11) {fksym=J I' fo^^ds 



-t/2 



belong to S^I^_^{T* X) for \t\ < TEhr, and this time is sharp. 

Let us insist on the fact that controlling the time evolution up to times t x 
log 1//?. is a crucial ingredient in order to obtain refined spectral estimates on Anosov 
manifolds. This will also be the case in Sections 13.31 and (3.41 when proving hyperbolic 
dispersion estimates and spectral gaps. 

The pseudodifferential calculus on '^1^^_^{X) allows to extend the validity of the 
bound (13. 6p up to the Ehrenfest time. The large deviation estimate (13.91) then leads 
to the following fractal Weyl upper bound for the number of atypical quantum decay 
rates. 

Theorem 14. [AnanlOj ylss'Mme the geodesic flow on S*X is Anosov. Then, for any 
a > q and e > 0, one has for h small enough 



(3.12) # <! 2 G Spec P{h) n D{l/2, Ch), >a \ <h-i 



-#^ + (l-rf)-€ 



where H is the rate function ( [J. 7[ j. A similar expression holds when counting quan- 
tum decay rates smaller than a' < q. 

Since —H{a) > for all a > q, this upper bound improves the bound o{h~'^^^) of 
Thm. [TT]by a fractional power of h (the name "fractal Weyl upper bound" takes its 
origin in the counting of resonances of chaotic scattering systems [SjoZwoOTp . For 
any a G [q, q+] one has H{a) > — sup^/i(v9+) > — Amax('^ — 1)? so the exponent of h 
in (I3.12P is negative, allowing the presence of (many) quantum decay rates arbitrary 
close to g+. 

These fractal upper bounds are not expected to be sharp in general [AnanlOj . In 
the next sections we will be able, under certain conditions, to exclude the possibility 
of quantum decay rates near q^ (or near g_). 



^For any fc G E, (5 g [0, 1/2), the symbol class Sg{T*X) consists of functions g{x,(,; h) which 
may become more and more singular when fi — >■ 0, but in a controlled way: 



In this class one can still use pseudodifferential calculus, and the expansions in powers of h make 
sense |EvZwi Sec. 4.3]. 
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3.3. A pressure criterium for a spectral gap. The result of Thm [9l namely 
the fact that all quantum decay rates belong to the interval [g_ — o(l),g+ + o(l)], 
was obtained by studying the norm of the propagator = e~**^/'' for times |t| ^ 1 
independent of h. Indeed, the sequence of equahties fl2.5p directly leads to the upper 
bound 

(3.13) Im Zn/h < q+ + e. 

The crucial semiclassical ingredient in fl2.5p is the second equality connecting the 
norm of the PDO B{t)x{Po) with the supremum of its (principal) symbol. Hence, 
one cannot hope to improve the bound fl3.13p as long as B(t)x{Po) remains a "good" 
PDO, which is the case if the time t < TEhrf^- To improve on f l3.13p we will 
investigate the propagator V^* for "large logarithmic times" , namely 

(3.14) t ~ /Clogl/^, with /C > large (but independent of h). 

For such times, the function b{t)x{po) oscillates on scales much smaller than h, so 
it cannot belong to any decent symbol class. As a result, the norm of B(t)x{Po) is 
a priori unrelated with that function. However, using the hyperbolicity of the flow 
one can prove an upper bound for the propagator in terms of a certain topolog- 
ical pressure depending on the damping and the hyperbolicity, and then obtain a 
constraint on the quantum decay rates. 

Before stating our bound, let us first introduce the notion of topological pressure 
associated with the flow $* on S*X. For any observable / G C{S*X), the pressure 
V(f) = p(/, \s,x) can be defined by 

(3.15) V{f) = sup {hKsif^) + /^(/), f^eM}. 

Notice that the definition (13. 7p of the rate function H{s) is very similar. The pressure 
is also given by the growth rate of weighted sums over long closed geodesic^: 

P(/) = lim log V e-^y^. 

T<h\<T+l 

The pressure appearing in the next theorem is associated with the function / = 
q — /2, where (p^ is the infinitesimal unstable Jacobian (13. 3p . it thus mixes the 
damping and hyperbolicity. 

Theorem 15. [SchelOj -Le^ X be an Anosov manifold, and q G S^(T*X) a damping 
function. 

Assume the following (purely classical) inequality holds: 

(3.16) P(g-V?+/2)<g+. 

Then, for any e > 0, there exists N > such that, for h small enough, the following 

resolvent bounds: 

(3.17) 

G D{l/2,Ch) n {Imz/h > V{q - (/?+/2) + e} , \\{P{h) - zy^\\ < fT^ . 



^The sum over long closed geodesies is analogous to a partition function in statistical mechanics, 
with the time T corresponding to the volume of the system. This explains why the growth rate is 
called a "pressure" . 
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Figure 3.2. Spectral gap in case the pressure V{q — '{'^ < q+. 
As a consequence, all quantum decay rates for ZniK) G D{l/2,Ch) satisfy 

(3.18) lHi|M<p(,_^./2) + . 

in particular we have a spectral gap. 

In the next subsection we present situations for which the inequahty (13.161) is 
satisfied. The proof of this theorem was inspired by a similar result in the case of 
resonances in chaotic scattering [NZ2[ INZ3] . We give some hints of the proof in 

Let us now apply this bound to the damped wave equation, that is take q{x,^) = 
—a{x). 

Corollary 16. |SchelO] Consider the DWE on X with nontrivial damping a > 0. 
Assume that the GCC fails (a- = 0), but that the conditionV{~a — (f^ /2) < holds 
true. 

Then the DWE has a spectral gap, and the energy decays exponentially for regular 
initial data, as in U.10\) . 

3.3.1. Conditions for a spectral gap. The inequality f l3.16p can hold only if the damp- 
ing function q varies "sufficiently" around its average value q. Indeed, one always 
has V{—ip^/2) > for an Anosov fiow, so by adding a constant dampint q = q one 
gets 

Viq - ^+/2) = q + P(-<^+/2) > q = q+. 

The pressure depends continuously of q (in the topology), so if q = q + 6q with 
||5g||^o small, the inequality (I3.16P will not be satisfied. 

Let us consider the case of the DWE as in Corollary [161 The failure of the GCC 
is equivalent to the fact that the set of undamped trajectories K (11.91) is nonempty. 
This set is fiow-invariant and closed, so one can define the pressure V{—(p~^/2, $* \x) 
associated with the fiow $* restricted on K. It was shown in [Schell] that if we 
multiply the damping a{x) by a large constant C > (this does not modify the set 
K), then 

lim V{-Ca - ip+/2) = V{-ip+/2, $* \k). 

C— )-oo 

If the set K is thin enough, the pressure V{—(f'^/2, $* \k) is negative (for instance, if 
K consists in a single closed geodesic, the pressure is equal to the average of —ip~^/2 
along the geodesic, which is negative). In that case, the above limit shows that the 
pressure V{—Ca — ip^ /2) is also negative ii C ^ 1. 
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3.3.2. Proof of the pressure bound: decomposing into "symbolic paths". In this sec- 
tion we explain the strategy of proof of Thm [151 As mentioned before, we want to 
bound the norm of the propagator for large logarithmic times. We will prove the 
following 

Lemma 17. Choose e > 0. Then, for an energy cutoff xo of small enough support, 
and a large enough constant /C > 0, for h < one has the bound 

(3.19) ||V^*Xo(i'o)|| < e*(^(^-'^^/'^+^), t~/Clogl/n. 

From there, obtaining the resolvent estimate fl3.17p is relatively straightforward. 

Proof. Our aim is to bound the norm || V^^mqH for a normalized state Uo microlocalized 
in Pq"^ ( [1/2 ± 5] ) , and a time t ^ IC log 1/h. We may assume that Uq is microlocalized 
in an open subset W of diameter < 5. We will first decompose mq in an adapted 
basis of local "momentum" modes {cr,, rj G M*^, \ri\ < 5}. Here the parameter rj is 
the "momentum coordinate" in a Darboux coordinate frame {(i/,?]), \y\ < 5, |?7| < 5} 
in ly, chosen such that r^i = — 1/2 represents the energy, yi the time, and the 
"momentum" Lagrangian leaves A^;, = {(y, ?7o)} are close the unstable foliation. The 
momentum states Cr, are Lagrangian states supported by the leaves A^, of norms 
0{1). The decomposition of Uq in this "basis" reads 

(3.20) Uo = Mv) e, j^:^ + o{hn, 

where the factor (2iTh)~'^^'^ ensures that HmoIIli = C^(l)- The reason for this de- 
composition is that we are able to precisely control the evolution of each individ- 
ual state Cfj, up to times t ~ /Clogl//i. The state y*e^ is a Lagrangian state 
supported on the Lagrangian leaf $*(A^). Due to the hyperbolicity of the flow, 
when t grows the leaf ^'(A^) expands exponentially fast along the unstable direc- 
tions, and converges to the unstable foliation. To analyze the state V^ejj, we can 
view as the combination of many local Lagrangian states associated with bounded 
pieces of ^^(A^). A convenient bookkeeping consists in using a partition of unity 
in Pq^ ([1/2 ±5]). One first covers this region by a finite family of open subsets 
{Wj, j = 1, . . . , J} of diameters < 6, and then constructs a smooth partition of 
unity {vTj G C^(Wj, [0, 1]), j = 1, . . . , J} adapted to this cover: J2j=i — ^ ^^^^ 

Pq^ ([1/2 ± S]). This partition is quantized into jllj =^ Opu^nj), j = 1, . . . , j|, such 
that, for any energy cutoff x any supported in [1/2 ± 5], 

J 

(3.21) J]n,x(Po) =x(^o) + 0L2^L2(^'"). 

We use this quantum partition to split the propagator V"" into "symbolic paths" : 
J 

(3.22) y%= Va,...a„er, + 0{h°^), V^,...^^ = U^^VU^„_,---VU^,V. 

ai,...,a„=l 

For each symbolic sequence ct = ai - ■ ■ an, the state V^e^ is a Lagrangian state sup- 
ported on the Lagrangian leaf A" obtained from by a sequence of n evolutions by 
$^ and truncations on Wa^- The full evolved leaf is given by $*(A^) = IJ|a|=nA"- 
If A" is empty V^e^ = 0{h°°). Otherwise {a. "admissible"). A" is a leaf of di- 
ameter < 6, close to the unstable foliation. One can compute the amplitude of 
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the Lagrangian state V^e,, to any order in h. It is governed by both the accumu- 
lated damping and accumulated instability along the "path" a, leading to the norm 
estimate 

(3.23) \\Vc,eJ<Cb{a)J+{a)-^^\ where 

n 

(3.24) 6(a) = TT ^(afc), 6(j) = max6(l, p), and similarly for (j^) ^^'^ . 

fe=i 

Applying the triangular inequality to (I3.22p . we obtain a bound on the norm of l^"e^ 
which can be related with the topological pressure: 

J 

(3.25) ||V%||<C J2 Kf^)J^{ay^^^<ew{n{V{q-^+/2) + 0{6))}. 

a admis. 

Inserting this bound in the expansion fl3.20p we get the same bound for V^uo with an 
extra factor h~'^^'^; however, this factor is smaller than e"*^ if we take n ~ ATlogl/Zi 
and K, large enough, in which case the bound (13.251) also applies to ||\^"uo||. This 
ends the proof of the norm bound (13.191) . □ 

Had we inserted (I3.23P in the expansion (I3.20p before summing over a, we would 
have obtained the following hyperbolic dispersion estimate: 

IIK.II <min(6(a), C h-^^^Hict) J+{a)-^/^) , \a\ </Clogl/n. 

This type of estimate was first proved in the case of the undamped operator Po{h) 
on an Anosov manifold |Anan08l lANlj . The adaptation to the case of the damped 
wave equation was written in [SchelOj . A similar estimate was also shown in [NZ2j 
in the case of chaotic scattering. 

3.4. A "thickness" condition for a spectral gap. In this section we present a 
spectral gap condition expressed only in terms of the set of undamped trajectories 
(II. 9p . but independently of the strength of the variations of the damping (as op- 
posed to the "pressure condition" (13.161) ). A proof of this result will appear in a 
forthcoming publication. 

To state our result in the context of a damping function q{x,^), we introduce the 
set of "least damped trajectories" 



(3.26) K = [J {snppn, IJ, e M, iJ,{q) = q+}, 

which generalize^ the set (11.91) . Our condition for a spectral gap depends on a 
topological pressure associated with the flow $* \k- We use quantities defined in 

mi 

Theorem 18. Let X be an Anosov manifold, and q G S'^{T*X) a damping function. 
Assume the following condition holds: 

(3.27) P(-y.+ , \k) <{d-l){^- A,nax) . 

Then, there exists 7 > ("the gap"), C > and N > such that, for h > small 
enough, the following resolvent bound holds: 

(3.28) yzeD{l/2,Cfi)n{lmz/h>q+--f}, \\{P{h) - zy^\\ < Ch,-^ . 



'''The two definitions are not strictly equal [Schellj . but this subtlety will be irrelevant here. 
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As a consequence, all eigenvalues of P{h) in D{l/2,Ch) satisfy 

Imzrjh < g+ - 7, 

so we have a spectral gap. 

In the framework of the damped wave equation with (nontrivial) damping function 
a(x) > 0, this gives the following 

Corollary 19. Assume that the set of undamped trajectories U.y\) is such that the 
condition \3.21\) holds. Then, there exists 7 > 0, C > and N > such that 

\\{T-A)-'\\^,^^,<Cr'', VrG{|Rer| >C, Imr>-7}. 

As a consequence, the energy decays exponentially for regular initial data, as in 

On a manifold of constant negative curvature — A^, the condition fl3.27p takes the 
form 

(3.29) htopi^' \k) <{d- l)A/2, 

where htop is the topological entropy of the flow on K. It is usually interpreted as a 
measure of the "complexity" of the flow on K\ we rather see it as a measure of the 
"thickness" of K. This bound reminds us of the lower bound obtained by Anan- 
tharaman when describing the localization of eigenstates of the Laplacian on Anosov 
manifolds [AnanOSj : in constant curvature, she shows the semiclassical measures as- 
sociated with any sequence of eigenstates cannot be supported on a set of entropy 
smaller than [d — l)A/2, thus forbidding the eigenstates from being too localized. 
Following the thread of Remark [T21 the spectral gap condition fl3.29p could be seen 
as a "nonperturbative analogue" of this delocalization result. 

Remark 20. For certain Anosov manifolds of dimension (i > 3, the condition (13.271) 
can never be satisfied. Indeed, the variational formula fl3.15p shows that, for any 
closed invariant set K' C S*X, the pressure 

\k') > - sup /i(^+) > -{d - l)A^ax. 

iJ.eM 

If the instability along is very anistropic (meaning that the largest and smallest 
positive Lyapunov exponents are very different), then {d — l)Amax — ^'^P^i^m 
can be strictly larger than i^til^^jsnijn ^ which case 

V{—(p^, \k') > {d — 1) (^—^ — Amaxj for any closed invariant set K'. 

This remark hints at the fact that the condition fl3.27p is probably not sharp 
(even in 2 dimensions), except maybe on manifolds of constant curvature. Riviere 
has recently improved Anantharaman's lower bound on the support of semiclassical 
measures in variable curvature |Rivll] : he shows that any such support S must 
satisfy V{—ip~^/2, $* \s) > 0. If we follow the above analogy (and also following the 
results of §3.3.ip . it seems natural to expect the following 

Conjecture 21. // the set K of least damped trajectories satisfies 

(3.30) P(-(/.+/2,$* t^) <0, 

then there is a spectral gap in the spectrum of P{h), in the sense of ThmlTEl 
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In constant curvature the condition fl3.30p is equivalent with f l3.29p . In variable 
curvature, it is weaker than fl3.27p . A Proof of f l3.30p should make use of local 
expansion rates, instead of the globally defined rates Amax and Umm, like in Riviere's 
work on 2-dimensional Anosov manifolds |RivlO] . 

In the next subsection we sketch the proof of Thm [181 

3.4.1. Sketch of proof for ThmUR As was the case for Thm[T5| the proof proceeds by 
bounding the norm of the propagator V^XoiPo) fo^' some logarithmic time (now the 
double of the Ehrenfest time, see fl3.38p ). In the course of the proof we will also need 
to control the evolution of a certain type of "microscopic" Lagrangian states. The 
novelty compared to the proof of Thm[T5]is that we will now distinguish between two 
types of phase space points, the "weakly damped" vs. "strongly damped" points. 

Without loss of generalit}!!, we may assume that the damping function q{x,^) is 
compactly supported inside p^^ ([1/2 ± e]). From now on we denote the Ehrenfest 

time by T =^ Tshr ■ We fix some level a G {q,q+), and consider the set of "weakly 
damped points" 

^+,a =^ {P G T*X, {q)T,sym{q) > «} 

(remember that {q)t,sym is the symmetric time average (13. lip ). The large deviation 
estimate (13.91) provides a bound on the volume of Q+^a H S*X: 

(3.31) fiL{^+,a n s*x) < /I- w-^(^). 

From the definition (13. 7p . H{q+) is equal to the pressure P(— $* \k) appearing 
in (I3.27p . If we assume (I3.27p . then, by continuity of H{s) on g+j, we may choose 
a G {q,q+) large enough such that 

(3.32) /3{a) = {d - 1) f ^ - J - ^(a) > . 



We will now associate a quantum projector to the set f2+,a- We first symmetrize the 

factorization (12.30 . by writing 

(3.33) 

V* = f/*/25^(t)f/</2^ where Bs{t) has principal symbol ^^(t) = e*^^^*'^'"". 

Although this symbol is positive, the operator Bgit) may not be selfadjoint; we then 
take its polar decomposition 



Bsit) = W{t)A{t) , where 



A{t) = {Bs{tYBs{t)f''^ is definite positive, 
W{t) is unitary. 



For the same reasons as in ^3.2.31 (and due to the support assumption on g), 
the operators Bs{T), A{T), W{T) remain "good" PDO up to the Ehrenfest time: 
Bs{T), A{T) G n-^vI/?/2_e(^), and W{T) G ^/^.^X). 

The operators A{T) and Bg(T) have the same leading symbol bs{T) = e-^^'^^^-^?'™. 
Hence, to the set fl+^a = {bs{T,p) > e"^} we associate the spectral projector 

n+ = n+.Q, =^ ]lyi(T)>eQT , and call n_ = I — 11+. 



'Lemma [7] can be extended to logarithmic times. 
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We use these projectors to decompose the propagator at time 2T: for some energy 
cutoff xo we write 

(3.34) V'^Xo(Po) = U^'^W{T) A{T) U^W{T) A{T) U^'\o{P^) 

= u^/^w{T) A{T) (n+ + n_) u^w{T) (n+ + n_) a{t) u^^\o{Po). 

The RHS sphts into four terms. Three terms contain at least one factor n_: for 
them we use the obvious bound 

(3.35) ||A(T)n_|| = ||n_A(T)|| < e"^. 
The remaining term contains the factor 

ul^ = u+u'^w{T) n+. 

The norm of this operator will be bounded by the following hyperbolic dispersion 
estimate, the proof of which is sketched in §3.4.21 

Proposition 22. Assume that for some a G (g, q+) the condition i\3. 32^) holds. Fix 
some e > 0. Then, if the energy cutoff Xi has small enough support, for h> small 
enough one has 

(3.36) \\Ul^XiiPo)\\,:^,.<h^-''^'\ 

Since propagation does not modify the energy localization, if choose in fl3.34p a 
cutoff Xo -< Xi) we then have 

A{T) U^'\o{Po) = XiiPo) A{T) U^I\,{P,) + 0(/i°^). 

Inserting this identity and the bounds f l2.7p . f l3.35p and f l3.36p in the identity fl3.34p . 

we get 

(3.37) \\V^^Xo{Po)\\ < e2^(^++^(^)) (e'^^^") + e^^^-''^) + e^^^^-^^)) . 

We may optimize this upper bound over the level a: the optimal value of the 
exponent is reached for the (unique) parameter ac G (g, g+) solving 

/3(q;) = g+ — a. 

For any 7 > satisfying g+ — 7 > we get (for > small enough) the 

following norm bound for the propagator: 

(3.38) ||^^^Xo(Po)||<e^^(^+-''). 

The proof of the resolvent estimate fl3.28p is then rather straightforward. 

3.4.2. Proof of the norm hound for f/J+. In this last subsection we sketch the proof 
of Proposition [22l that is obtain an upper bound for 

(3.39) \\uI_,xi{Po)\\l2^l^ = ||n+f/^iy(T)n+xi(i^o)||, T = TEhr. 

We recall that 11+ is the projector associated with the region i7+,Q,. Let us indicate 
that a similar type of dispersion estimate was used by S. Brooks, when studying 
the delocalization of the eigenstates of quantized hyperbolic automorphisms of the 
2-dimensional torus (the so-called "quantum cat maps") [BrolOj . 

To estimate this norm, it will be useful to replace this projector by a smoothed 
microlocal projector obtained by quantizing a symbol x+ = X+,a ^ S^^^_^{T*X), 
such that Opji{x+) "dominates" 11+: 

n+xi(^o) = Op,(x+) n+xi(^o) + 0(^°°)- 
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The norm fl3.39p can then be bounded by: 

\\Ul^Xi{Po)\\ < sup \{0p,ix+)u,,U^WiT)0p,{x+)ui)\+O{h^). 

Il"l|| = ll«2||=l 

The symbol x+ can be chosen supported inside a set of the form Q_cel~'pQ ^([l/2± 
e]). This set is quite irregular, and the main information we have on it is an estimate 
on its volume (using large deviation estimates like fl3.3ip ). It is then convenient to 
use an anti-Wick quantization scheme for Opfc(x+), that is use a family of coherent 
states (Gaussian wavepackets) {cp, p G T*X}|j, to define 

dp 

Each coherent state Cp is normalized, and is microlocalized in a "microscopic ellipse" 
around p. This ellipse is chosen to be "adapted" to the flow. Let us describe it 
using the local Darboux coordinates {{y,Ti)} as in §3.3.21 The ellipse is "short" in 
the energy direction, At]i ~ h}~'^^'^, and "long" in the time direction, Ayi ~ h^^"^. 
The spread along the transverse directions is chosen isotropic: Ayj = Ar]j ~ h}^'^, 
j = 2,...,d. 

The scalar product (Op^(x+)M2, ?7^VFs(T) Op'^{x+)'^i) can now be expressed as 
a double phase space integral 




dpi dp2 
Po'([i/2±e]) (2vrn) 



2d («2,ep,) {e,„Ui)x+{P2)x+{Pi){U~^^'e,„U^/^W{T)e,,). 



The state W(T)ep-^ is approximately identical to Cp^. We can then precisely describe 
the (undamped) evolutions of the coherent states Cp- up to the times ±T/2. The 
state tZ-^/^Cpj (resp. U~'^^'^ep^) is a (microscopic) Lagrangian state along a leaf of 
the weak unstable manifold of volume ~ h^'^^^^'^^/'^ J^j^ipi) centered at $"^/^(pi), 
(resp. a leaf of the weak stable manifold of volume ~ h'^'^^^^'^^/'^ J^i^{^~'^/'^{p2)) 

centered at $~"^/^(p2))- Using the sharp energy localization of these states and the 
fact that stable and unstable manifolds intersect transversely to each other, one gets 
the following bound: 

Q f m {pi)-m{p2) \ 

^J+,(pi)4/2($-^/2(p2)) 

for some 9 G C^([— 1, 1]). The denominator is bounded below by e'^'^in{d-i)T/2 ^ 
Inserting this bound in the above double integral and using large deviation estimates 
similar with (13.311) for the support of x+, one finally gets (after some manipulations) 
the bound (IX^ . □ 
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